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RESEARCH MEMORANDUM 

LONGITUDINAL STABILITY, TRIM, AND DRAG CHARACTERISTICS 
OF A ROCKET-PROPELLED MODEL OF AN AIRPLANE 
CONFIGURATION HAVING A 45° SWEPTBACK 
WING AND AN UNSWEPT HORIZONTAL TAIL 
By James H. Parks and Alan B. Kehlet 


SUMMARY 


A free-f light investigation was made of the longitudinal stability, 
trim, and drag of a rocket-propelled airplane model at low lift coef- 
ficients at a range of Mach numbers from O .63 to l.l 6 . The configu- 
ration included a wing and horizontal tail of aspect ratio 4 and 
thickness ratio in the streamwise direction of 6 percent. The quarter- 
chord line of the wing was sweptback 45°, whereas the tail was unswept 
and fixed at 2° incidence, trailing edge down. The center of gravity 
was located at zero percent of the mean aerodynamic chord. Oscillations 
induced by pulse rockets were used to obtain stability data. 

The lift-curve slopes were linear over the range tested and reached 
a maximum value of O.O 96 at M = 0.91- Calculations based on structural 
influence coefficients obtained from static tests indicate appreciable 
losses in lifting ability especially at the higher Mach numbers. The 
configuration exhibited a high degree of stability and stable damping 
characteristics over the speed range investigated. 

A smooth nose-up trim change of low magnitude occurs near M = 0.90. 
The drag coefficient at trim lift increases from a subsonic value of 
0.021 to 0.055 near M = 1.10. 

INTRODUCTION 


A general research program has been initiated by the National 
Advisory Committee for Aeronautics to determine, by means of rocket- 
propelled models in free flight, the effect of various empennage designs 
on the transonic longitudinal stability, trim, and drag characteristics 
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of complete airplane configurations. Presented herein are the results 
from one of the transonic models having a conventional empennage 
arrangement with an unswept horizontal tail fixed at 2° incidence, 
trailing edge down, and mounted on a 45° sweptback vertical tail. The 
wing was sweptback 45° with an aspect ratio of 4 and thickness ratio 
of 6 percent in the streamwise direction. 

Longitudinal stability, trim, and drag were obtained from an 
analysis of continuous telemeter records and of short-period oscilla- 
tions induced by vertically thrusting pulse rockets. 

The model was tested at the Pilotless Aircraft Research Station 
at Wallops Island, Va. 


SYMBOLS 


exposed span, ft 


c 

g 

P 

q 

y 

K 

% 

t 

L 

M 

P 

R 

S 


mean aerodynamic chord, ft 
gravitational acceleration, ft/ sec 2 
static pressure, lb/sq ft 
dynamic pressure, "L pM^, lb/sq ft 

increment of exposed span, ft 

ratio of lift-curve slope of a rigid wing to that 
of present wing 

radius of gyration about y-axis, ft 

time, sec 
unit load, lb 
Mach number 

period of the short-period oscillations, sec 
Reynolds number, based on c 
wing area, sq ft 


CONFIDENTIAL 


NACA RM L52F05 


CONFIDENTIAL 


3 
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W 

I y 

t 1/2 

a n/s 

a Z/6 

a 


60 

7 

9 

C N 

C L 

C D 

Cm 

c m 0 

Subscripts : 
a 


q 


velocity, ft/sec 
weight of model, lb 

moment of inertia about the y-axis, slug-ft 2 
time to damp to one -half amplitude, sec 
normal accelerometer reading, positive up 
longitudinal accelerometer reading, positive forward 
angle of attack, deg 

horizontal tail deflection for zero angle of attack 
ratio of specific heats, 1.4 
angle of pitching, radians 

angle of twist, radians; leading edge positive up 
normal force coefficient, (w/s) (l/q) (a n /g) 
chord force coefficient, (w/s) ( l/q) ( -a i/g) 
lift coefficient, cos a - C c sin a 

drag coefficient, -C c cos a - sin a 
pitching-moment coefficient 

pitching-moment coefficient at zero angle of attack 


(da/dt)(c/2V)(l/57.3) 

(de/dt) (c/2V) 


The symbols a, a, and q used as subscripts indicate the 
derivative of the quantity with respect to the subscript; for example, 
C L a = dC L /da. 
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MODEL AND APPARATUS 
Model 


A three-view drawing of the model is shown in figure 1(a). The 
model was constructed mainly from laminated mahogany. Metal plates 
incorporated in the wing and horizontal tail for additional stiffness 
and rigidity are shown by the sectional details of figure 1 (b). 

The wing had an aspect ratio of 4, taper ratio of 0.60, and 
NACA 65 AOO 6 airfoil sections in the streamwise direction with the 
quarter-chord line sweptback 45°. The horizontal tail had the same 
geometrical characteristics except the quarter-chord line was unswept. 

The vertical tail had an aspect ratio of 1.5, taper ratio of 0.50, and 
NACA 65 AOO 8 airfoil sections in the streamwise direction with the 
quarter-chord line sweptback 45°. The wing had no incidence, whereas 
the horizontal tail was fixed at 2 ° incidence, trailing edge down. 

The center of gravity of the model was located at zero percent of the 
mean aerodynamic chord. 

The fuselage was a parabolic body of revolution of fineness ratio 
8.91 which is described fully in reference 1. Fuselage ordinates are 
given in table I. 

Photographs of the model are shown in figure 2. A small metal 
hook of l/ 8 -inch steel was attached to the fuselage for boosting 
purposes. The relative size and location of this hook is shown clearly 
in figure 2(b). The leading edge of this hook had a sharp wedge section. 


Propulsion 


The model-booster combination is shown on the launching platform 
in figure 3- The model was launched at an angle of 60° elevation and 
was boosted to maximum velocity by an ABL Deacon rocket motor. 


Six vertically thrusting pulse rockets were installed in the model 
(four in the nose section and two in the rearward section) . The 
locations of the pulse rockets are shown in figure 1(c). Each pulse 
rocket had a total impulse of approximately 8 pound- seconds and a 
burning time of approximately 0.08 second. 
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Instrumentation 

The model was equipped with an NACA four-channel telemeter which 
transmitted continuous records of normal and longitudinal accelerations, 
angle of attack, and total pressure. 

The flight path was determined from tracking radar data and 
atmospheric conditions at altitude were obtained from a radiosonde 
released immediately after model firing. 


DATA ANALYSIS 


The technique of data reduction for an analysis of the response 
of models to abrupt disturbances is described in reference 2 for abrupt 
elevator deflections. The method applies equally well for models 
employing pulse rockets. Briefly, however, static longitudinal stability 
is determined from the periods of the short-period oscillations and 
dynamic longitudinal stability is determined from the rate of decay of 
the oscillations. The oscillations occurring during pulse rocket 
burning are not included in the analysis because the time history of 
the thrust-forcing .function cannot be evaluated accurately. 

The trim lift, angle of attack, and drag were determined between 
pulses directly from the telemetered data and through oscillations by 
appropriate fairing. The angle-of -attack data were converted to angle 
of attack of the center of gravity by the methods of reference 3 • 

A detailed discussion of the accuracy of this type of investigation 
is found in reference 2. For the particular instrumentation used, the 
absolute accuracy of and C-q is ±0.010 and ±0.003* respectively, 

at M = 1.10, and ±0.025 and ±0.008, respectively, at M = 0.70. The 
angle of attack is believed to be correct within 0.20° and the Mach 
number is estimated to be correct within 0.02 at M = 1.00. 

Since the influence of aeroelasticity on aerodynamic characteristics 
is important when sweepback is incorporated in a lifting surface, 
static tests were made to determine the flexibility of the exposed 
portion of the model wing. Concentrated loads were applied independ- 
ently at several spanwise stations along the 25 -percent and 50 -percent 
chord lines and the resulting angles of twist along the span for each 
loading condition were determined. Deflection diagrams obtained from 
these data are presented in figure 4. The variation of dynamic pressure 
with Mach number is shown in figure 5 for use with structural influence 
coefficients in the calculation of aeroelastic effects. 
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RESULTS AND DISCUSSION 


The Reynolds numbers of the test are shown in figure 6 as a 
function of the Mach number. A portion of the flight time-history 
converted to aerodynamic parameters is shown in figure 7. Seven 
oscillations, similar to those shown in figure were available for 
analysis; six induced by pulse rockets between Mach numbers of 0.92 
and 0.65 and one at separation (M = 1.16) resulting from the dif- 
ferences in trim between the model alone and the model-booster com- 
bination. The last pulse rocket fired before the previous oscillation 
had completed a sufficient number of cycles for the usual stability 
analysis and thus is included in the lift analysis only. It should be 
pointed out that all the data analyzed were in the lift-coefficient 
range of ±0.20 as shown in figure 8. 


Lift 

The lift curves obtained are plotted in figure 8. The values of 
lift-curve slope represented by the faired lines in figure 8 are shown 
in figure 9- The subsonic value of 0.079 increases abruptly to 0.096 
near M = 0.90. Near a Mach number of 1.11, the lift-curve slope is 
0.06l. These values are less than would be obtained from a similar 
configuration having a rigid wing because of the effects of flexibility 
mentioned previously. 

By using the deflection data of figure A in the form of influence 
coefficients, as suggested in reference h, a factor by which the 
flexible wing data may be corrected to the rigid wing case was computed. 
This factor is shown as a function of Mach number in figure 10. It 
should be pointed out that this correction is not precise because of 
the type of data available but rather is presented to indicate the 
order of magnitude of the losses due to flexibility. 

Values of lift-curve slope corrected for the effects of wing 
flexibility are shown in figure 11 j also shown are lift-curve slopes 
calculated from the wind-tunnel tests of component parts reported in 
references 5 and 6. The agreement is good except in the region between 
M = 0.90 and M = 1.00. The reasons for the discrepancy in this 
Mach number range are not completely known but similar effects have been 
noted previously between free-flight and wind-tunnel tests. 


CONFIDENTIAL 


MCA RM L52F05 


CONFIDENTIAL 


7 


Static Longitudinal Stability 

Periods of the short-period oscillations are shown as a function 
of Mach number in figure 12. Since no pulse data were obtained between 
M = 0.91 and M = 1.11, a dashed-line fairing is used over this 
relatively large Mach number increment in the stability analysis plots. 
These dashed-line fairings are based on unpublished rocket-model data 
and general considerations of the present test. These data converted 
to the static stability parameter Cn^ are shown in figure 13 . The 

configuration is shown to be longitudinally stable throughout the Mach 
number range for the center-of -gravity location used. In general, the 
longitudinal stability increased with increasing Mach number. Although 
some loss occurred between M = 0.80 and M = 0.87, no severe varia- 
tions were noted and the loss had been regained at M = 0.91* 

The degree of longitudinal stability, as indicated by the 
aerodynamic -center location, is shown in figure 14. At Mach numbers 
less than 0.87, the shape of the curve is quite similar to the 
C n ^ variation. 


No attempt was made to isolate the factors contributing to these 
variations as the influence of wing flexibility was also felt at the 
tail in the form of changes in downwash. Generally speaking, however, 
the two major effects of wing flexibility on the configuration stability 
are in opposite directions; the loss of lift over the wing combined 
with a forward shift in the wing center of pressure tend to decrease 
the over-all stability, whereas the moment -producing ability of the 
tail should increase due to changes in downwash, thereby tending to 
increase the stability. 


Dynamic Longitudinal Stability 

The times required for the short-period oscillations to damp to 
one-half amplitude are shown in figure 15* Damping factors corresponding 
to these time increments are shown in figure 1 6 . The damping is stable 
over the Mach number range tested. 

Although the total damping factor decreases from approximately 60 
at high subsonic Mach numbers to 48 at M = 1.11, the damping-moment 
derivative remains essentially the same, approximately 25. No attempt 
was made to isolate the variables contributing to the damping caused 
by the effects of wing flexibility, but apparently the moment-producing 
ability of the tail increased, whereas the over-all lifting ability, 
supplied mainly by the wing, decreased with increasing Mach number. 
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It should be pointed out that a mass-distribution term Ky/c was 

included in the total damping factor. For the present tests, 

Ky/c = I. 89 . The importance of this term in comparing the present 

results with the two-degree-of -freedom damping characteristics of con- 
figurations with different mass distributions is discussed fully in 
reference 7- 


Longitudinal Trim 

The trim lift coefficients for the configuration through the Mach 
number range are shown in figure 17(a) • The shape of the trim curve 
agrees with the results of reference 8. The fact remains that a pulse 
rocket fired during the transonic trim change may obscure the exact 
nature of the variation, but the transition from subsonic to supersonic 
flight appears to introduce a smooth nose-up change in trim with the 
trim lift coefficient increasing from 0 at M = O .85 to O.Ok at 
M = 0.95. 

Angles of attack corresponding to these trim lift coefficients are 
shown in figure 17(b). Although the changes in lift-curve slope are 
reflected in these data, no unusual variations were noted. 

By reasonably assuming C m to vary linearly with Cp (from basic 

oscillations), a trim curve for a center -of -gravity location at c/k 
was computed and is shown in figure 17(a)* The decrease in stability 
amplifies the magnitude of the trim change but the shape of the trim 
curve remained essentially the same. 

On a similar basis, the trim angles of attack were used to 
determine the pitching-moment coefficient at zero angle of attack. 

These values are shown in figure 18. As might be expected from the 
data presented previously, these pitching-moment coefficients increased 
with increasing Mach number with the most abrupt change in the region 
of M = 0.90. 

It is interesting to note that if the horizontal tail is considered 
as an all-movable control, zero angle of attack could be maintained 
throughout the Mach number range with little control movement as shown 
in figure 19 . These values were computed from the relationship 
&o = C mo /C m5 where C m& is from the wind-tunnel data of reference 5* 

The variations in C mQ and 5 0 are believed to be independent 

of the wing and primarily due to Mach number effects on the down flow 
over the horizontal tail induced by the convergence of the rearward 
part of the fuselage from the subsequent considerations. The 
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wing-fuselage combination is symmetrical about the longitudinal axis 
and no downwash from the wing should exist at zero angle of attack. 

The contribution of the empennage drag to the zero-lift pitching-moment 
coefficient is calculated to be relatively small, 0.003 at M = 0.70 
and 0.012 at M = 1.10. The same effects are shown in the more complete 
data of reference 9 wherein similar total changes of approximately 0.06 
in CniQ and 1° in horizontal tail deflections were noted. 


Drag 

The drag coefficients at trim lift coefficients throughout the 
test Mach number range are shown in figure 20. The subsonic value of 

0.021 increasing to 0.055 at M = 1.11 with the most abrupt increase 
near M = 0.95 is of the order of magnitude which might be expected 
from the geometry of the configuration. The reasons for the hook and 
subsequent increase in abruptness of the drag rise near M = 0.9^ is 
believed due to pressure changes on the rear fuselage and have been 
noted previously on similar configurations (ref. 10). 

Minimum drag values are also shown in figure 20 and, except near 
M = 1.11, the trim lift had no appreciable effect on the drag coef- 
ficient. The lift-coefficient range of the transient oscillations was 
not great enough to determine the effect of lift on drag with any degree 
of accuracy. 


CONCLUSIONS 


From the flight tests at low lift coefficients of a rocket- 
propelled model of an airplane configuration having a 45° sweptback 
wing and an unswept horizontal tail with 2® incidence, the following 
conclusions may be drawn: 

1. The lift-curve slopes were linear over the ranges tested with 
a maximum value of 0.096 at M = 0.91- Calculations based on static 
loading data indicated that losses in lift due to wing llexibility 
occurred over the entire speed range with the magnitude of the losses 
increasing with increasing Mach number. 

2. The aerodynamic -center location remained near 60 percent c 
at subsonic Mach numbers. At M = 1.11, however, the aerodynamic 
center location had moved rearward to 88 percent c. 

3. The configuration exhibited stable damping characteristics 
over the Mach number range. 
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t. A longitudinal transonic trim change occurred near M = 0.90 
in a nose-up direction. The change in trim was smooth and of small 
magnitude. 

5. The drag coefficient at trim lift increased from 0.021 at 
subsonic speeds to 0.055 at supersonic speeds. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I 

FUSELAGE ORDINATES 


Station 
(in. ) 

Diameter 

(in. ) 

0 

0 

3 

1.60 

6 

3.00 

9 

4.24 

12 

5.28 

15 

6.14 

18 

6.84 

21 

7-34 

24 

7.66 

27.8 

7.80 

30 

7.78 

33 

7.74 

36 

7.64 

39 

7.48 

42 

7.30 

45 

7.06 

48 

6.78 

51 

6.44 

54 

6.08 

57 

5.66 

6o 

5.18 

63 

4.68 

66 

4.12 

69.5 

3.42 
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(a) General arrangement of model. 


Figure I.- Arrangement and construction details of model. 

are in inches. 
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(b) Details of wing and tail cpnstruction. 


VJ1 


Figure 1.- Continued. 
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(c) Arrangement of the pulse rockets in the fuselage. 
Figure 1.- Concluded. 
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(a) Top view. L"72331*l 



(b) Side view. 

Figure 2 .- Photographs of the model. 



L-72332.1 
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(b) Loads at 0.50 chord line. 

Figure 4.- Deflection diagrams obtained from static loading tests on the 

exposed model wing. 
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.6 .7 .8 .9 1.0 1.1 1.2 

Mach number , m 

Figure 5.- Variation of dynamic pressure with Mach number. 
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Figure 6.- Variation of Reynolds number with Mach number. 
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Figure 7*- Portion of the time history of the flight converted to 

aerodynamic coefficients. 
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Mach number, M 


1.16 - 1.08 


.92 - -90 


.89 - .86 
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.65 - .63 


- 4 - 2 0 2 4 

Angle of attack, a, deg 

Figure 8.- Lift-curve data obtained from the short period oscillations. 
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Figure 9.- Variation of the lift-curve slope with Mach number. 
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Figure 10.- Factors for converting the flexible wing 

to the rigid case. 
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Figure 11.- Comparison of lift-curve slopes from two techniques with the 

effects of aeroelasticity considered. 







j 



















O'" 







° 1 
















— 

- Pres 
Wine 

sent t 
1 turn 

:es ts 
lel t e 

x K 
ss ts 

( Hefs , 

• 5 ai 

id 6) 














i 

.NACA 




Figure 12.- Variation of the periods of the transient oscillations with 

Mach number. 
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Figure 13*- Variation of the static stability parameter with Mach number. 
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Figure 14.- Variation of the aerodynamic-center location with Mach number. 
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Mach number , M 
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gure 15.- Times required for the short period oscillations to damp to 
one -half amplitude as a function of Mach number. 



Mach number , m 


Figure 1 6 .- Variation of damping factors with Mach number. 


CONFIDENTIAL 




Trim angle of attack, deg Trim lift coefficient 


28 


CONFIDENTIAL 


NACA RM L52F05 



.6 .7 .8 .9 1.0 1.1 1.2 

Mach number , M 


(a) Lift coefficient. 



(b) Angle of attack. 

Figure 17.- Longitudinal trim characteristics as a function of Mach number. 
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Figure 18. - Variation of the pitching-moment coefficient at zero angle of 

attack with Mach number. 
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Figure 19.- Horizontal tail deflection required to maintain zero angle 
of attack ^determined by using C m& from ref. 5^. 



Figure 20.- Variation of drag coefficients with Mach number. 
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